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i ^ ■ Abstract. The time delay effect for planets and spacecraft is obtained from a fully 

relativistic modified gravity theory including a fifth force skew symmetric field by fitting 
to the Pioneer 10/11 anomalous acceleration data. A possible detection of the predicted 
CN| ■ time delay corrections to general relativity for the outer planets and future spacecraft 

^ \ missions is considered. The time delay correction to GR predicted by the modified 

gravity is consistent with the observational limit of the Doppler tracking measurement 
reported by the Cassini spacecraft on its way to Saturn, and the correction increases 
I to a value that could be measured for a spacecraft approaching Neptune and Pluto 

o 

o 

O" 1 . 1. Introduction 

U 

In previous work [1], we have provided a possible gravitational explanation for the 
Pioneer 10/11 anomalous acceleration, based on the scalar-tensor-vector gravity (STVG) 
theory [2]. In the following we will obtain predictions for the corrections to the GR 
(general relativity) Shapiro time delay [3, 4, 5, 6, 7], obtained from the fitted anomalous 
acceleration observed in the Pioneer 10/11 spacecraft data [8, 9, 10, 11, 12]. For our 
study of the solar system, we must account in our modified gravity theory for the 
variation of G with respect to the radial distance r from the center of the Sun. Since 
we do not possess rigorous solutions of the field equations, we use a phenomenological 
parameterization of the varying parameters a(r) and A(r) to obtain fits to the anomalous 
Pioneer acceleration data that are consistent with the solar system and fifth force 
experimental bounds. Thus, with the variation of G with distance from the Sun, we 
can make predictions for other observational tests in the solar system. In the following, 
we shall calculate the corrections to the GR time delay effect for planets and spacecraft 
probes obtained from the STVG modified gravity theory. 

2. Pioneer Anomalous Acceleration 



We assume that gravity is the cause of the Pioneer 10/11 anomaly due to the difference 
between the running G(r) and the bare value, G?o ~ G^, where G^ denotes the 
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Newtonian gravitational constant. The Pioneer anomalous acceleration directed towards 
the center of the Sun is given by [1, 2]: 

5G{r)M & 



ap 



where 



SG(r) = G N a{r) 1 - exp(-r/A(r)) ( 1 + 



A(r) 



(1) 



(2) 



We will use the following parametric representations of the "running" of a(r) and 
A(r) [1]: 

a(r) = a 0O (l-exp(-r/f)) fc / 2 , (3) 

A(r) = A 00 (l-exp(-r/f))- fe . (4) 

Here, f is a non-running distance scale parameter and b is a constant^. 

A best fit to the acceleration data extracted from Ref. [11] has been obtained using 
a nonlinear least-squares fitting routine including estimated errors from the Doppler 
shift observations [9]. The best fit parameters are [1]: 

aoo = (1.00 ±0.02) x 10~ 3 , 
Aoo = 47 ± 1 AU, 
f = 4.6±0.2AU, 

b = 4.0. (5) 

The small uncertainties in the best fit parameters are due to the remarkably low variance 
of residuals corresponding to a reduced x 2 per degree of freedom of 0.42 signalling a 
good fit. 



3. Time Delay in Modified Gravity 

Our scalar-tensor-vector gravity (STVG) theory does not have an exact spherically 
symmetric static solution. However, we have shown that the solution will have the 
form of the Schwarzschild solution for large values of the radial coordinate r [2]. In the 
derivation of the line element, we have neglected the "running" of G. The line element 
in isotropic coordinates is given by 

(-> gm\ 2 

ds 2 = ) J - 2 c 2 dt 2 - (l + — ) da\ (6) 

/ 1 , GM \ v ic rJ 

\ L 2c 2 r J 

where 

da 2 = dr 2 + r 2 d9 2 + r 2 sin 2 Odtf. (7) 
% Note that in ref. [1] the exponent b in Equation (12) should be —b. 
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Expanding in powers of GM/r, we obtain in Cartesian coordinates 

2GM /GM\ 2 i „ (_ 2GM* 



ds z 



c 2 r 



(-5-) c 2 dt 2 - 1 + — 5— )(dx 2 + dy 2 + dz 2 ), (8) 
V c r / J V c r / 



where r = (x 2 + y 2 + z 2 ) 1 / 2 ,^ = arctan[z/(x 2 + y 2 ) 1 ^ 2 ] and <p — arctan(y/x). The Sun 
is taken to be at the origin of coordinates and Earth and the spacecraft (planet) lie in 
the z = plane and the transmission null ray lies along the x direction. For a null ray 

ds 2 = g 00 c 2 dt 2 + g u dx 2 = 0. (9) 

We now replace G in (8) by G^ + 5G where 5G is given by (2) and obtain the correction 
to the GR round-trip time delay: 



St = / dx 



5G((x 2 + rg) 1 /*)-, 



(x 2 + r 2 )V2 

where r ffi and r p are the distances of the Sun from Earth and the spacecraft (planet), 
respectively, and ro = y = const. The correction 5r in (10) is obtained by performing a 
numerical integration. The excess time delay for a round-trip delay is given in GR by 

(Ar) GR = —^\n(—). (11) 

The excess delay is calculated for a signal that grazes the limb of the Sun, r ~ Rq, and 
it is a maximum when the spacecraft (planet) is at superior conjunction. 

In ref. [1], the variation of 5G/Gn arising from Equation (2) for the parametric 
values of a(r) and A(r) of Equation (3) and Equation (4), respectively, was obtained 
using the best fit values for the parameters given in Equation (5). The behavior of 
G{t)/Gn is closely constrained to unity over the inner planets until beyond the orbit of 
Saturn (r > 10 AU) where the deviation of G(r) from Newton's constant increases to 
an asymptotic value of Goo/G^ — > 1.001 over a distance of hundreds of AU. 

In Table 1, we display the predicted values of the GR time delay (At)gr, the 
correction 5t obtained from STVG and the observational limit. 



Table 1. Theoretical predictions for the corrections St to the GR time delays 



compared to the observational limits. 


r 


GR prediction STVG prediction 


Observational limit 


(AU) 


(At) gr (pisec) 


ST(fisec) 


( At) obs limit (/isec) 


1.52 (Mars) 


247 


3.07 x 10~ 12 


0.5 


5.20 (Jupiter) 


271 


3.10 x 10~ 7 




8.43 (Cassini) 


262 


8.71 x 10~ 6 


6.0 x 10~ 3 


19.22 (Uranus) 


297 


3.49 x 10~ 4 




30.06 (Neptune) 


305 


1.15 x 10~ 3 




39.52 (Pluto) 


311 


2.06 x 10~ 3 




100 


329 


9.32 x to" 3 





The Cassini spacecraft while on its way to Saturn reported the most accurate 
Doppler tracking measurement of the time delay effect [6]. The accuracy of the 
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measurement was made possible by using both X-band (7175 MHz) and Ka-band (34316 
MHz) radar, thereby reducing significantly the dispersive effects of the solar corona. 
The 2002 superior conjunction of Cassini was favorable with the spacecraft at 8.43 AU 
from the Sun, and the distance of closest approach of the radar signals to the Sun was 
r ~ 1.6 R Q . 

We see from the result for the correction to the GR time delay for the Cassini space 
probe that the prediction is well below the observational limit of 6.0 x 10~ 3 /zsec. The 
reported value for the measured post-Newtonian parameter 7 is 7 — 1 = (2.1±2.3) x 10~ 5 . 
Forthcoming new data for the anomalous acceleration will provide a more accurate 
determination of the parameters for a(r) and A(r) and the predicted values of St. 

A spacecraft orbiting Neptune and Pluto or on its way out of the solar system could 
hopefully report a sufficiently accurate Doppler tracking measurement of the time delay 
that could detect the correction to the GR prediction. At the positions of Neptune and 
Pluto and for a spacecraft at a distance from the Sun, r = 100 AU, the corrections 5r 
to the GR time delay are of the order of the observational limit of the Cassini Doppler 
tracking measurement, so Doppler measurements with comparable accuracy reported 
by a spacecraft approaching Neptune and Pluto and beyond the solar system could 
distinguish GR from the modified gravity theory prediction based on the Pioneer 10/11 
data. 

The Doppler measurements should be carried out at conjunction with the spacecraft 
(planet) to maximize the effects of the curvature of space at the distance of closest 
approach to the Sun. Moreover, Doppler measurements should be reported by the 
spacecraft (planet) at opposition to determine a baseline that yields the Newtonian 
time taken by a signal to reach the spacecraft (planet) and return to Earth. This latter 
timing will require using the JPL ephemeris including the modified gravity (STVG) to 
determine the corrected orbital parameters of the spacecraft (planet). Subtracting the 
baseline timing from the measurements at the distance of closest approach to the Sun 
will result in a residual correction St that can test the modified gravity theory. It should 
be noted that the calculation of 5t is not significantly sensitive to increasing values of 
r - 

4. Conclusions 

We have demonstrated in ref. [1] that the STVG theory can explain the Pioneer 
anomalous acceleration data and still be consistent with the accurate equivalence 
principle, lunar laser ranging and satellite data for the inner solar system as well as the 
outer solar system planets including Pluto at a distance of r = 39.52 AU = 5.91 x 10 12 
meters. The ephemerides for the outer planets are not as well know as the ones for the 
inner planets due to their large distances from the Sun. The orbital data for Pluto only 
correspond to the planet having gone round 1/3 of its orbit. It is important that the 
distance range parameter lies in the region 47 AU < A(r) < 00 for the best fit to the 
Pioneer acceleration data, for the range in the modified Yukawa correction to Newtonian 
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gravity lies in a distance range beyond Pluto. 

Perhaps, a future deep space probe can produce data that can check the predictions 
for the corrections to the time delays obtained from the Pioneer anomaly acceleration 
data and from our modified gravity theory. An analysis of anomalous acceleration data 
obtained from earlier Doppler shift data retrieval will clarify in better detail the apparent 
onset of the anomalous acceleration beyond the position of Saturn's orbit. 
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